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ABSTRACT 



Experimental results for the convective heat transfer and flow friction 
characteristics of plate-fin type heat transfer surfaces are presented for 
eight surfaces. Six surfaces were fabricated of perforated nickel plate, 
one perforated nickel fins with solid nickel plate splitters and one of 
solid nickel plate. The heat transfer data were obtained by the transient 
technique and includes the effect of longitudinal conduction. 
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NOMENCLATURE 



English Letter Symbols 



A 


Matrix total heat transfer area, ft^ 


A 

c 


2 

Matrix minimum free flow area, ft 


A. 

fr 


2 

Matrix total frontal area, ft 


A 

s 


Solid matrix cross sectional area available for thermal 
conduction, ft^ 


\ 


2 

Conduction area corrected for effect of perforations, ft 


A* 


2 

Plane surface area, ft 


a 


Plate thickness, ft 


a 


Short side of a rectangular flow passage, ft 


b 


Plate spacing, ft or in. 


b 


Long side of a rectangular flow passage, ft 


C 


Flow stream capacity rate, (mc^) , Btu/fhr 8 F) 


C 

s 


Matrix capacity, W g c g , Btu/°F 


c 

P 


Specific heat at constant pressure, Btu/(lbm°F) 


c 

s 


Matrix material specific heat, Btu/(lbm°F) 


d h 


Hydraulic diameter, 4r^ 


d 


Inside pipe diameter, in. 


E 


2 

Flow friction power per unit area, HP/ft 


f 


Mean friction factor, dimensionless. This is the "small” 
or "Fanning" friction factor. (Ratio of wall shear stress 
to the fluid dynamic head.) 


G 


* 2 

Exchanger flow stream mass velocity, (m/A c > , lbm/(hr ft ) 


8 C 


Proportionality factor in Newton's Second Law, g = 32.2 
(lbm ft)/(lbf sec^) 


h 


Unit conductance for thermal convection heat transfer, 
Btu/ (hr ft z °F) , or heat transfer power per unit area 
per degree temperature difference, Btu/(hr ft^ °F) 
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English Letter Symbols (continued) 



j 


2/3 

Colburn factor * Ng T , heat transfer characteristic, 

dimensionless 


j' 


Colburn j factor based on plane surface area, A*, dimension- 
less 


K , K 
c ’ e 


Contraction loss coefficient for flow at heat exchanger 
entrance or exit respectively, dimensionless 


k 


Unit thermal conductivity, Btu/(hr ft^ °F/ft) 


k 

s 


2 

Matrix thermal conductivity, Btu/ (hr ft °F/ft) 


L 


Total matrix flow length, ft 


% 

m 


Mass flow rate, lbm/hr 


P 


2 

Pressure, lbf/ft 


P 


Porosity for matrix surfaces, dimensionless 


q 


Heat transfer rate, Btu/hr 


R 


Gas constant, (ft lbf)/(lbm °R) , (53.35 for air) 


r h 


Hydraulic radius, (A^L/A) , ft, (4r^ * hydraulic diameter) 


s 


Solidity of a perforated material, 1-u, dimensionless 


T 


Absolute temperature, degrees Rankine, °R 


t 


Temperature, degrees Fahrenheit, °F 


V 

m 


3 

Matrix volume, ft 


V 

s 


Material volume corrected for effects of perforations 


W 

s 


Mass of matrix, lbm 


X 


Distance along the flow passage in direction of flow, ft 



Greek Letter Symbols 



a 


Aspect ratio of a rectangular flow passage, b/a, dimension- 
less 


6 


Compactness; ratio of total heat transfer area to the volume 
between the plates, ft 2 /ft^ 


A 


Deno tes difference 


6 


Fin thickness, ft 
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Greek Letter Symbols (continued) 



C 

0 

X 

T 

M 

u 

p 

0 

♦ 

1 

Z 

Subscripts 

atm 

f 

1 
k 
m 

0 
s 

STD 

1 

2 



Conduction area reduction ratio due to perforations in 
perforated material; cross sectional area of perforations/ 
cross sectional area, dimensionless 



Time 

Longitudinal conduction parameter, 



k A 



s s 
ALc 



p 



dimensionless 



Time parameter 

Fluid viscosity coefficient, lbm/hr ft 

Area reduction ratio due to perforations in perforated 
material, dimensionless 

3 

Density, lbm/ft 

Ratio of freeflow area to frontal area, A /A, , dimensionles 

c fr 

Denotes "function of" 



Compactness for perforated material including the effect 

of area reduction, A*s/V , dimensionless 

m 

Summation 



Local atmosphere 
Fluid 

Initial, individual 
Equivalent 

Mean or matrix, as appropriate 
At orifice 

(Matrix material) solid 

Standard temperature and pressure 

Upstream 

Downstream 
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Dimensionless Groupings 



N Reynolds number, (4r,G/y) , a flow modulus 

K n 



N. 



R(p) 

N, 



St 



N 



Nu 



Pr 



N, 



Tu 



Reynolds number for pipe, (dG/y) 

Stanton number, (h/Gc^) , a heat transfer modulus 

Nusselt number, (h4r^/k) , a heat transfer modulus 

Prandtl number, (yc^/k) , a fluid properties modulus 

Number of heat transfer units, (hA/mc ) 

P 

Generalized heat transfer grouping - Colburn "j" factor, 
(Nst^pr ' /• This factor versus defines heat transfer 
characteristics of the surface. 



f 



X 



X 



k 



T 



Mean friction factor. This is the "small 11 or "Fanning" 
friction factor (Ratio of wall shear stress to the fluid 
dynamic head). This factor versus Nr defines the fric- 
tion characteristics of the surface. 

Conduction parameter, k g A s /mLc p for solid material; 
kgAk/mLcp for perforated material 

Equivalent conduction parameter (^orrected for equivalent 

length in perforated material) , Xj^ 

k 

h A 

Time parameter, ~ 0 

s s 
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INTRODUCTION 



The present Investigation is an outgrowth of the continuing program 
at Stanford University [8]^ for the evaluation of heat transfer surfaces. 
The experimental work reported by Howard [A] was performed at the U. S. 
Naval Postgraduate School, Monterey, and a continuation of that work is 
presented herein. 

DESCRIPTION OF SURFACES 

The eight compact heat exchanger surfaces tested were: 

1. Skewed passage 

a. 20° total skew angle; 160/A0 TV perforated nickel 

2. Modified parallel plate 

a. 160/40 TV perforated nickel 

b. 160/40 Q perforated nickel 

c. 125 M perforated nickel 

d. 125 P perforated nickel 

e. 50 G perforated nickel 

f. 160/40 TV perforated nickel fins with solid nickel splitters 

g. solid nickel 

The skewed passage surface employs corrugated sheets of 160/40 TV per- 
forated nickel, stacked alternately with the lines of corrugations 10° 
from the flow direction, thereby forming a total angle of 20° between 



1 



Numbers in brackets refer to numbered items listed in References. 
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adjacent stacked plates. It is the same geometry as one employed by 
Howard [4]. The modified parallel plate surfaces were made by alternately 
stacking a formed plate and a plane splitter plate to prevent nesting, as 
shown in the sketch below and in the enlarged photograph of Figure 30. 



The flow channels formed by this geometry may be compared to rectangular 

channels with an aspect ratio, a = b/a, of about 7. 

The perforated nickel plate employed is an electro-deposited metallic 

sheet of integral structure manufactured by Perforated Products, Incorporated. 

This material is described in detail in Appendix C. The geometrical and 

physical properties of the surfaces are given in Figures 3 through 7 . 

All matrices formed had a frontal cross section of approximately 

3.2 inches square and a length of 2.0 inches in the flow direction. All 

matrices had flow passages with a hydraulic diameter of approximately 

.002 feet, thus the L/D u ratio in all cases was of the order of 83. 

ri 
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PRESENTATION OF RESULTS 



Heat transfer and flow friction data for each matrix investigated 
are given in both tabular and graphical form, TABLES I through VIII and 
Figures 9 through 16, using the Colburn j modulus. Fanning friction factor, 
and Reynolds number. In evaluating the Reynolds number, the hydraulic 
diameter, which is defined as four times the hydraulic radius, was used. 

The effects of entrance, exit, and flow acceleration have been considered 
in the evaluation of f. (See Appendix A) . 

Figures 17 through 22 show heat transfer and friction comparisons 
between various matrices investigated. Also plotted in Figure 18 for 
comparison purposes are the theoretical laminar flow solutions for 
parallel plates and rectangular passages [6, 8, 11]. Both solutions are 
for an infinite length to hydraulic diameter ratio, and the heat transfer 
solutions are for constant wall temperature. In Figures 23 through 25, 
"figure of merit" type curves are presented. The j/f versus Reynolds 
number thus plotted gives an indication of the required matrix flow 
frontal area for a given pressure drop. 

Lastly, in Figures 26 through 28, heat transfer power versus flow 
friction power curves on a unit area basis, evaluated for fluid proper- 
ties at standard conditions of dry air at 500°F and one atmosphere, are 
presented (See Appendix A). 

EXPERIMENTAL METHODS 

To obtain the heat transfer data, a transient technique was employed. 
Briefly, this method consists of heating the test matrix by heated air 
to a uniform temperature (approximately 20°F above ambient in these tests) 
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and then subjecting the matrix to a step change in the air flow tempera- 
ture to ambient temperature , The air temperature downstream of the 
matrix is monitored and recorded versus time. There are several ways 
of constructing such a test rig [4, 10, 12] • The one used in these tests 
is shown in Figure 1, and described in detail in Appendix B. 

A series of fine nichrome wire heaters were installed upstream of 
the matrix to heat the air, Turning off the heaters provides the step 
change in temperature of the air without disturbing the flow. By 
referencing the thermocouples downstream of the matrix to the thermo- 
couples upstream of the heaters, the initial temperature difference can 
be very closely controlled. This temperature difference is continuously 
recorded by a strip chart recorder when a run is made. This recorded 
trace has the distinct advantage in that there is no transposition of 
data required which would produce increased uncertainties. Reduction of 
the data follows the method of Locke [10] and Howard [5], whereby the 
of the surface can be evaluated by determining the maximum slope of 
the fluid temperature difference versus time curve during the cooling 
transient. 

Pressure drop data for evaluating the friction factor were obtained 
from static pressure taps located in the test section immediately upstream 
and downstream of the matrix. The static pressure immediately upstream 
of the matrix was also recorded. Velocity profiles were taken just 
upstream of the matrix to assure that uniform flow conditions were being 
maintained by the screen pack straightener and by the wire heaters. The 
flow was measured by an A 0 S.M,E, standard D and D/2 orifice meter with 
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changeable orifice plates [14], and was located downstream of the test 
section. Various manometers were used consistent with the pressure 
range encountered. All pressure drop data for evaluating the friction 
factor were taken under ambient air conditions (i.e., isothermal flow). 
The data reduction relationships are given in detail in Appendix A. 

EXPERIMENTAL UNCERTAINTIES 

The experimental uncertainties can be considered to be in two cate- 
gories. First, in the manner in which the experimental gear attempts to 
meet the idealizations of the experimental method. These are discussed 
in detail by Howard in [4] and [5]. 

By insuring uniform flow to the matrix, minimizing the temperature 
change to minimize fluid property changes, and in the careful design 
of the fluid heater, these idealizations are met. The second catagory 
can be grouped in three areas : 

(1) uncertainty of physical constants, 

(2) inaccuracy in the geometric measurements , and 

(3) instrumentation inaccuracies. 

For the determination of the uncertainties reported, the method of Kline 
and McClintock [9] was used. 

(1) Values for the physical constants required were obtained from 
[1], [2], and [3]. The uncertainties in these values, as best as can be 
determined, are listed below: 

c * 0.5% 

s 

c * 0.5% 

P 

N_ * 2.0% 

Pr 

k * 0.5% 

s 

p * 1.0% 
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It is noteworthy that the accuracy for the value of C g depends upon 
the metal used* Where pure metals are used the accuracy quoted is 
feasible; however, this is not the case with alloys* Thus the tolerance 
of c g can be a major cause for inaccuracy [10]. 

(2) Due to inconsistencies in construction of the matrices, errors 
in lineal dimensions are considered to be ± 0*5%* Inasmuch as the 
weight of the matrix can be determined as accurately as desired, the 
error in W g is considered negligible* Based on this, the uncertainties 
in geometric measurements are considered to be as follows: 

A, ^c * ^s * ’ ± 

L ± 0o 5% 

W < 0*1%, negligible 

s 

(3) Instrumentation inaccuracies were essentially those obtained 
in pressure measurement* With the exception of the temperature measure- 
ment at the orifice, all temperatures recorded were temperature dif- 
ferences which were recorded in inches (See Appendix A)* The temperature 
measurement at the orifice was obtained by a copper-constantan shielded 
thermocouple which was read by a Leeds and Northrup portable potentio- 
meter. Assuming adequate manufacturers calibration of the thermocouple 
wire, the estimated possible error is ± one-half of the smallest division 
on the potentiometer, or ± 0*0025 millivolts which is approximately ± 0.1°F* 

Since the range of pressures varied, requiring several different 
manometers, the maximum uncertainty encountered was used in the uncertainty 
analysis. These values are listed below: 

P ± 1**25% 

o 

AP ± 1 o 25% 

o 

AP ± 1*70 % 

m 

P atm 4 0o0005 " H S (negligible) 
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One final parameter for which uncertainty must be determined before 
an analysis of the experimental results can be made is N,^. Inasmuch as 
Colburn j factor is not a linear function of Reynolds number, uncer- 
tainty for high and low N values was considered. Using Figure 3-A 
of reference [4] or Figure 5 of reference [5] and the uncertainty in 
maximum slope as 2.0% the uncertainty interval using reference [9] 
in the Colburn j factor varies from * 7.5% at high Reynolds Numbers to 
* 10.2% at low Reynolds Numbers. The uncertainty interval for Reynolds 
Number is * 2.3%, and for the friction factor, f, the uncertainty 
interval is * 4.3%. 

DISCUSSION OF RESULTS 

The perforated nickel matrices investigated covered a wide range of 
percent open area, 12 to 50%. Both slotted and round perforated nickel 
surfaces were evaluated for the "modified parallel plate" geometry. The 
perforated nickel geometric properties are given in APPENDIX C for the 
various types investigated. In order to make a reasonable comparison of 
the performance among these matrices, the same frontal area and flow 
length were specified. Results are presented in Figures 9 and 10 and 
listed in TABLES I and II for the slotted perforated nickel matrices; 
those for the round perforated matrices are given in Figures 11, 12, and 
13 and tabulated in TABLES III, IV, and V. 

In order to establish a meaningful fiducial point, a solid nickel 
matrix of the identical material (obtained from the same manufacturer) 
with the same frontal area and flow length as the perforated matrices was 
fabricated. The performance of this matrix is illustrated in Figure 15 
with tabular values given in TABLE VII. 
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The best performing slotted perforated nickel matrix was the 160/40Q. 
The 50G round perforated nickel (50% open area) was the best surface in- 
vestigated among the n parallel plate" nickel matrices. 

A$ another approach to evaluating the Colburn j factor, j’ was con- 
sidered. The parameter j f is based on the plane surface area, A*, as 
if no perforations were present; whereas, j is based on the heat transfer 
surface area. A, which includes the plate solidity correction, s. The 
value of j f is readily attainable from the relationship j* - j x s, 
where s = A/A*. (See Appendix D). The determination of the solidity 
is fully described in APPENDIX C. Figure 13 for the 50G perforated 
nickel matrix includes the values for j 1 versus Reynold’s number in 
addition to the f and j values. 

A comparison of the performance of the two slotted perforated 

% 

matrices together with the solid nickel matrix is shown in Figure 19. 

A similar comparison for the round perforated nickel matrices is given 
in Figure 20. 

It is readily apparent that the perforated material yields a 
higher heat transfer characteristic. Inasmuch as the friction factor 
also shows an increase, a more appropriate comparison can be made from 
the flow area "goodness" factor curve in Figures 23 through 25. The 
heat transfer power (hg TI) ) versus flow friction power ( E STD ) per unit 
area given in Figures 26 through 28, is another "goodness" factor 
evaluation. Here the information for each of the matrices is presented 
based on standard conditions for fluid properties and a common hydraulic 
diameter. (See APPENDIX A). 

So that an appreciation for the perforated material might be realized 
in direct transfer type of heat exchanger application in addition to its 
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use in matrix type for rotary regenerators, a "parallel plate" matrix 
with the 160/40TV perforated nickel fins and solid nickel splitter plates 
was fabricated*. The performance of this matrix is given in graphical 
and tabular form in Figure 14 and TABLE VI. A comparison with the solid 
nickel matrix is made in Figure 21 . Figure 25 contains the j/f presenta- 
tion and Figure 26 shows the vs. E gxD ^ or t * 1 * s ma trix. Here again 

is a noticeable increase in heat transfer performance, but with a 
negligible increase in friction factor* 

For the 160/40TV perforated nickel material, a 20° skew matrix was 
compared in performance against the "modified parallel plate" to see 
how these two configurations compare with one another. The separate 
performance for each has been mentioned previously; however, the 
direct comparison is shown in Figure 22* As with the other nickel surfaces, 
the flow area "goodness" factor is represented in Figure 23. 

In the results of this investigation are several items which are 
general in nature, but nonetheless extremely important. In the laminar 
Reynolds number range, friction factor is inversely proportional to 
Reynolds number. By Reynolds analogy, it was anticipated that the Colburn 
j factor should likewise be inversely proportional to Reynolds number in 
this region. However, the experimentally obtained values of j in the 
low laminar Reynolds range utilizing the equivalent conduction length 
(defined in APPENDIX C) together with Figure 2-A of [4] or Figure 4 of 
[5] yielded lower than anticipated values for j. Assuming the Reynolds 
analogy to hold the j vs. N experimental curve was extrapolated, shown 
as a dashed line on the Figures. The j and values obtained in this 

manner at the given experimental N are asterisked in the Tables. 

K 
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The second item is that of an equivalent conduction length which 
was stated earlier. In the perforated fin material the slots or round 
holes, as the case may be, provide for a winding heat flow path. It 
was appreciated that the matrix flow length, L, was not the true con- 
duction path but somewhat less than actual, A mean conduction path 
length evaluation is presented in APPENDIX C, The conduction parameter, 
X, based on the flow length, L, was evaluated, and by simply multiply- 
ing by the ratio L/L^ (also given in APPENDIX C) , the value of X^ was 
determined. 

The entrance and exit length effects for a matrix with a small 

hydraulic radius are small. In the evaluation for Fanning friction 

factor (APPENDIX A), the effects of exit and entrance length were 

considered using values of K and K from Figures 5-3, -4, or -5 of 

c e 

reference [8], as appropriate. For the perforated surfaces, the K values 

were taken corresponding to laminar flow conditions, 

A final item worth mentioning is again concerned with the low 

Reynolds range. A slight alteration in the conduction parameter value 

has a marked effect on the H value in this range. For this reason, 

when utilizing a nickel material, or any other material with a wide 

range of specified values for thermal conductivity, it is mandatory 

that the appropriate value of k g be obtained. The value of X is 

directly proportional to this k value, i,e,, 

s 



CONCLUSIONS 



A very definite improvement in performance is achieved with the 
utilization of a perforated material* Kays [6] postulated the 
improved heat transfer characteristic without "a large amount of 
form drag so characteristic of high performance surfaces". The 
perforations apparently disturb the thermal boundary layer to a 
much greater degree than the hydrodynamic boundary layer. The 
perforated nickel matrices investigated herein confirm Kays 1 
hypothesis. 

The experimentally determined Colburn j values in the low 
Reynolds range were lower than the j values obtained by extrapolating 
the j versus N curve from the higher Reynolds range. Because of the 
large uncertainties involved in the low Reynolds range the extra- 
polated data are considered more reliable. 

Effects of longitudinal conduction have been considered in 
the analysis of data. For the perforated material matrices, correc- 
tions for the actual conduction path length and for the area removed 
by the perforations have been included. 

The best performance for a perforated "parallel plate" matrix 
was attained by the 50G type material. 
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Figure 1. Schematic Diagram of Test Apparatus 
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Matrix Material 


160/40 TV 
Perf .Nickel 


160/40 Q 
Perf .NicKel 


Specific Heat (c s ) Btu/lb F 


0.1065 


0.1065 


Thermal Conductivity (k s ) 3tu/ nr ft F 


38.7 


38.7 


Material Thickness, inches 


0.0022 


0.0016 


Total Heat Transfer Area (A) ft^ 


18.593 


16.2667 


Frontal Area (Af r ) ft* 


0.06953 


0.06953 


Total Conduction Area (A k ) ft 2 


0.00424 


0.002887 


Free Flow Area (A c ) ft^ 


0.0584 


0.06142 


Matrix Vo lame ( t / Ti ) ft3 


0.011588 


0.011588 


Matrix Density ( ^ iTl ) lb/ft- 


63.6 


59.2 


Hydraulic Diameter (Dp) ft 


0.001929 


0.002028 


Compactness ( 0 ) ft*/ft3 


1604.5 


1403.8 


Porosity (p) 


0.840 


0.883 



3. Geometric *:.a Pnysical Properties of Slotted 
Perforated Mckei (160/40 TV ana I6O/4G 3) 
"Parallel Plate" Matrices 
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Matrix Material Perf . Nickel Perf. Nickel Perf . Nickel 
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Figure 4. Geometric and Physical Properties of Round Perforated 

Nickel ( 125 M, 125 P> and 50 G ) "Parallel Plate" Matrices 
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■ l 


Matrix Material 


160/40 TV Fins - 
Solid Nickel Splitters 


Specific Heat (c s ) Btu/lb °F 


0.1065 


Thermal Conductivity (k s ) Btu/' hr ft °F 


38.7 


Fin Thickness, inches 


0.0022 


Splitter Thickness, inches 


0.0020 


Total Heat Transfer Area (A) ft^ 


19.3084 


2 

Frontal Area (A^ r ) ft 


0.06953 


Total Conduction Area ( A^) ft* 


0.00689 


Free Flow Area (A c ) ft^ 


0.05883 


Matrix Volume (V m ) ft-' 


0.011588 


Matrix Density ( ^ m ) lb/ft^ 


69.1 


Hydraulic Diameter (D^) ft 


0.00194 


Compactness (Q) ft^/ft^ 


1666.2 


Porosity (p) 


0.846 



Figure 5. Geometric and Physical Properties of Slotted Perforated 
Nickel Fin (160/40 TV) with Solid Splitters,”Parallel 
Plate” Matrix 
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Matrix Material 


Solid Nickel 


Specific Heat (c s ) Btu/lb °F 


0.1065 


Thermal Conductivity (k s ) Btu/ hr ft °F 


38.7 


Material Thickness, inches 


0.0020 


Total Heat Transfer Area (A) ft 2 


20.1875 


Frontal Area (Af r ) ft 2 


0.06953 


Total Conduction Area (A g ) ft 2 


0.01014 


Free Flow Area (A c ) ft 2 


0.05939 


Matrix Volume (V m ) ft^ 


0.011588 


Matrix Density ( ^ m ) lb/ft ^ 


77.7 


Hydraulic Diameter (D^) ft 


0.001961 


Compactness ( Q ) ft 2 /ft^ 


1742.1 


Porosity (p) 


0.854 



Figure 6. Geometric and Physical Properties of a 
Solid Nickel "Parallel Plate 11 Matrix 




\ v A C \N- :. x V; V *v > >: 0 ^ x 

. Vs^'vy'' \ -V 0 AV\V ^) vi '"^Nv^v<;u 
■<■ v \ i y'£ ; r. i ; : j. ; j {■ 1 s > > * » '•* ; \ • > - » 

' • \ M, * ^ » It \ U v i i V t t » t i 



h VV? X. \\ \ 



> . , , , , i \ 

V v X,V v V. V v V, V 4 x A. v X . V x x v v x > * 
.' V "S *\ \ * • v\ vS k x.%vxx J V\ 



v * ' \ i \ ^ f V\ Cv X V ' ^ V ^ 4 ^ v ' ' x « \ X 

• H v\ ^\«\wAV\U * > }J >•£ >, > HA * > V A 

- - V 'A v * 'VV:-. V ^ lV ' Vl v vv v M * MV^^V\>VA \ A *i 

• v . x i f* » k ' ' > 41 x x v >.^*\ VV\:\V^ % 

} K * V \\ \ \\ X \ v * V V V * ? < '/ < * < X 



m.i v Os\\ • \ : v,v v y o- • >. 

\' s V V \ V \ > V y\ v V, ' *x \ X V V c X V x X \ v V 
v \ V y V 'O v m ^ > \ x t \ v v a \ v v v x \ \ \ 

', ' ' V V ' ' •>*<> |' ' v » x V 8 v v \ x s X A X X x x x X V 

X X , , ‘ ‘ 'x * 'X \ X X X X V \ x X x \ % A *\ A A ‘X X X X X x x X 

\ x -v\Vvr>.VV^V 0 \A& P \ \ * \k\ \ k 1 v > ' \ * ' * > > 

v v A A X \ A . vVx *v\ \\\ X N a v vA A \/x‘V*x A \ 



S' x s - v v v ^ x > X* w» V\ \ x \ x N A v Ox A \ A V\ \ \ 

\ V'> \ \\ V\ VY* * Y\Vv* V* i V&' . > VX . V.^'V 

$ . ' v n 5 * \ \ XWAX U 

V ' C '' V V x. v v' ' x x ' ' x' v \ X ' ' .4 x \ X Xr . x V < X X V x 

*■ x \ CS *x-^ ^^x ' V^* >'.- Vy>‘ V^ -'CSSi V ^ '*0 V , s , v V V •' X 

> s ^ > 'i > ''V^ > VV aA \ v\ w v.\ ?s\ XX \ \ \ \\ \ \ 

•v v « .VO W'v'^SstsSv O" v ' ' v " “ ‘A “ '> > 

t v O ‘ ' v , V \ * X \ X . V X x X x x l X X X « A > X x 



6 7 8 9 



THE L.B.STARRETT CO. Q 



1234'6789 



ATH0LMASS.U.S.A. 



Matrix Material 


20 Degree Skew 
160/40 TV 
Perf .Nickel 


Specific Heat (c s ) Btu/lb F 


0.1065 


Thermal Conductivity (k s ) Btu/ hr ft F 


38.7 


Material Thickness, inches 


0.0022 


Total Heat Transfer Area (a) ft^ 


15. a 


Frontal Area (A fr ) ft 2 


0.0734 


Total Conduction Area (A k ) ft 2 


0.00349 


free Flow Area (A c ) ft 2 


0.06412 


Matrix Volume (V m ) ft^ 


0.01211 


Matrix Density (^ m ) lb/ft 3 


50.1 


Hydraulic Diameter (D h ) ft 


0.00252 


Compactness ( ^ ) ft 2 /ft^ 


1272.8 


Porosity (p) 


0.873 



Figure 7. Geometric and Physical Properties of a Slotted 
Perforated Nickel (160/40 TV) 20° Skew Matrix 
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Figure 9. Surface Heat Transfer and Friction Data 

"Parallel Plate", 160/40 TV Perforated Nickel 




.Figure 10. Surface Heat Transfer and Friction Data 

"Parallel Plate", 160/40 Q Perforated Nickel 
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Figure 11, Surface Heat Transfer and Friction Data 

"Parallel Plate", 125 H Perforated Nickel 
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Figure 12. Surface Heat Transfer and Friction Data 

"Parallel Plate", 125 P Perforated Nickel 
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Figure 13. Surface Heat Transfer and Friction Data 
"Parallel Plate", 50 0 Perforated Nickel 
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Figure 14. Surface Heat Transfer and Friction Data 

"Parallel Plate", 160/40 TV Fins with Solid 
Nickel Splitter Plates 
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Modified Parallel Plate 
0.002 " Solid Nickel 



Porosity, p = 0.854 
Compactness, § = 1742.1 ft 2 /ft^ 
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Figure 15. Surface Heat Transfer and Friction Data 
"Parallel Plate", 0.002" Solid Nickel 
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Figure 16. Surface Heat Transfer and Friction Data 
20* Skew, 160/40 TV Perforated Nickel 








Figure 17. Comparison of Slotted Perforated Nickel, 160/40 TV vs. 
Solid Stainless Steel, 20° Skew Matrices 
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Figure 18. Comparison of Solid Nickel vs. Theoretical Solutions 
Parallel Plate Matrices 
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Figure 19. Comparison of Slotted Perforated vs. Solid 
Nickel "Parallel Plate" Matrices 







Figure 20. Comparison, of Round Perforated vs. Solid 
Nickel "Parallel Plate" Matrices 
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Figure 21. Comparison of Slotted Perforated Fin vdth Solid 

Nickel Splitters vs. Solid Nickel “Parallel Plate" 
Matrices 
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Figure 22. Comparison of Slotted Perforated Nickel "Parallel 
"i Plate" vs. Slotted Perforated Nickel 20° • Skew <• 

Matrix 
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Figure 23. Flow Area Goodness Factors for the 160/40 TV Perforated Nickel and the Solid 
Stainless Steel 20° Skew Matrices 
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Figure 27 » Heat Transfer Power vs. Flow Friction Power per Unit Area 
for the Round Perforated Nickel and Solid Nickel "Parallel 
Plate" Matrices 
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Figure 28. Heat Transfer Power vs. Flow Friction Power per Unit Area 
■ for v the 160/10 TV Perforated Nickel and the Solid 

Stainless Steel 20° Skew Matrices 



1 





\ 



■?> 



} 




Figure 29- Enlarged Photograph of the 20° Skew 
Matrix, 160/40 TV Perforated Nickel 




Figure 30. Enlarged Photograph of the Modified 
Parallel Plate Matrix, 160/40 TV 
Perforated Nickel 
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Figure 31* Matrix Holder 
(Upstream View) 
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© From Extrapolated Values of j versus No Curve 
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APPENDIX A 



DATA REDUCTION RELATIONSHIPS 



Summarized herein are those data reduction relationships which are 

of prime importance in obtaining the heat transfer and flow friction 

characteristics as presented from the raw data taken. Before the data 

on the matrices tested can be reduced, or a comparison of the different 

matrices can be made, it is essential to obtain accurate geometrical 

dimensions. The geometrical factors of interest are the porosity, p, 

the hydraulic diameter, D , and the area compactness, 6. Determination 

n 

of any of these fixes the third since it may be shown [11] that 



d h 

The definitions of these terms are: 

A 



(i) p = 



Free Flow Area 
Frontal Area 



fr 



(2) D„ = 



4 x Free Flow Area 
H Heat Transfer Surface Area 



= 4 -r- L 



(A-l) 

(A- 2) 
(A- 3) 



(3) S = 



Heat Transfer Surface Area 
Frontal Area x Flow Length 



A* L 
fr 



(A-4) 



MAXIMUM SLOPE 

The maximum slope of the generalized cooling curve is a unique 
function of and X [10]. 



t f2~ t i 



'tr 11 






- ») 



(A-5) 



max 



t and have both been previously defined, but for convenience 
the relationships are restated: 
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hA 



0 



and 



T 



W c 
s s 



Therefore : 



N, 



Tu 



hA 




t/n Tu - <v. c .> e - 

and 



d( T /N Tu ) = mc p /W g C g )d0 = (C/C g )d0 



(A-6) 



Furthermore : 



^2^1 

'n't 




fl 

i 



+ 



1 



and 



d 






i 




Combining (A-6) and (A-7) 



alM 

bfry 

d (t/n Tu i 



max 



t fr t i 



d[t f2~ t fl ] 

d0 



[7] 



(A-7) 



(A-8) 




Time — ► 
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From the above sketch 



(1) d[t f2~ t fl 1 



d0 



( 2 ) 



= X 
x 

max 

= d0, sec 



Chart speed 
(3) y = d(t f2 ~t fl ), inches 



(4) M = (t fl ~t i ), inches 



Combining these determined values into (A-8) yields 



4^1 

t fl~ t i i 
d[T/N Tu ] 



1 X 

— • -t— • 

M x 



Chart speed 



k A 
s s 

With this value of maximum slope and \ = v — — , enter TABLE I 

v me L 

P 

or Figure 4 of Howard [5] or Figure 2-A of Howard [4] to obtain a 

corresponding N value. 

tu 



(A-9) 



MASS RATE OF FLOW: 



The mass rate of flow is calculated from the method given in reference 

[14] with necessary rearrangements in order to reduce unnecessary hand 

calculations and to gain the full benefits of Murdock [13] : 

m = 359Kd 2 F Y/AP""y (A-10) 

o a o 

where: K = flow coefficient, including velocity of 

C 

approach = with C * orifice coefficient of 



discharge (obtained from [13]) and 

£^ = ratio of orifice diameter to internal pipe 

diameter, d Q /d 

d = orifice diameter, inches 
o 9 

= thermal expansion factor 
Y = expansion factor 

P 2 

y = specific weight of fluid flowing = — with P in lbf/ft . 

Ki 
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R = 53.35 (ft-lbf) / (lbm-°R) 

T = t(°F) + 459.7 = degrees Rankine (°R) 



AP q = Pressure drop across the orifice, inches 1^0 substituting 

p 

— for y, since the magnitude of the local acceleration of gravity, g, 

is taken as being equal to the magnitude of standard acceleration of 

c 

gravity, g ; and 1 ■ for K. gives: 
c rr7, 4 



% 

m 



359C ,2 
° 



F Yv4p • 
a o 



P 

RT 



(A-ll) 



Inserting the numerical value for R and converting P to units of 
lbf/in^ yields: 

: ° n " " 2 ” ” ' • £ (l>sla) (A-12) 



m = 



589.81 Cdp F a Y 



/ 1 - 6 4 






o T 



From Figure 38, reference [14]: F =: 1.00 

a 

From Figure 40b, reference [14]: Y = 0.9985 = 1.0 
The pressure, P, in (A-12), is hereafter referred to as P. 
and is determined by the following relationship: 



P, « (P - P /13.6) 0.4912, lbf/in 2 (A-13) 

3 atm o ’ 

Where P atm is the atmospheric pressure in inches of mercury, 

P q is the static pressure upstream of the orifice in inches 



of H 2 0 

13.6 inches of l^O/inch of mercury 
0.4912 (lbf /in 2 )/inch of mercury 

The temperature, T, in (A-12) is the air temperature at the 

orifice in degrees Rankine (T^) : 

T = t +459.7 (A-14) 

o o 

Replacing F and Y in (A-12) with their numerical values, the equation 
a 

reduces to the working form: 
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m = 



REYNOLDS NUMBER: 



589.81 C d 2 
o 



AP 



— lbm/hr 



(A-15) 



Reynolds number is defined as : 

d h g 
N = — — 

R p 



(A-16) 



where G is the mass flow velocity based on the free flow area, A^. 



G = 



m 



m 



pA 



c ‘ fr 

Substituting (A-17) into (A-16) : 

, .-i-i 

R > P 4S 1 

But from (A-l) D /p = — , therefore: N = — — — 

H p K ^ A f r ^ 



(A-17) 

(A->18) 
(A- 19) 



Equation (A-19) shows that for a given matrix and mass flow rate the 
Reynolds number is inversely proportional to the compactness. Substi- 



tuting (A-4) into (A-19) yields the working equation: 

N - 
R pA 



(A- 20) 



FANNING FRICTION FACTOR 



The following sketch together with equation (A-21) describes the 
flow system under consideration: [8] 

1 * b Z 



^ > 



(A-21) 

L.' ' ’ 0 • . 

v .1 iaoioq 



\ i \ * 

AP 
P, 



gLIi 

2 *c P 1 



(K +l-p 2 ) - 2 ( — -l)+f^- — - (l-p 2 -K ) — 
c,, W, J A v, e v- j 



entrance flow core exit 

effect acceleration friction effect 



' t 



i .17 tq- j. 



l/'j ' 



V / / / ( / / ////// 

i tezzzTZZzzzzzzh 
tzzzzzzzrzzzz zb 
tizzzzzzzzzzzzzh 

&ZZZZZZZZZZZZ2& 

fe zz zzzzzzzz zzh 
^zzzzzzzrzzzz^ \ 
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It Is important to note that v cs v. and v, ss v„ in the above 

a l t> z 

sketch as the pressure changes (for gas. flow heat exchanger application) 
from section 1 to a and b to 2, respectively, are very small relative 
to the total pressure. Inasmuch as the testing is performed with air 
at moderate temperatures and pressures, the perfect gas law is applicable 



(P = pRT) . 

^ D h A r fl Ac 

Substituting v = “ and recalling that — = — or — = — from 

equation (A-3) , equation (A-21) solved for f for the isothermal case 

becomes : 

f 



r, AP P l +P 2 f K C , K e\ P l +P 2/l 1_V 1+ 2\1 !h , a 
] 2g c pm q7 ~ 2 \P 1 + P 2 )” 2 (p 2 ~ P^ )\ 1+P /J L (A " 22) 



where pm 



p l +p 2 



with the upstream values subscripted by 1 and the 



downstream subscripted by 2. K c is the entrance coefficient and K g the 
exit coefficient, and both are dependent on porosity, shape of the flow 
cross section, and the matrix Reynolds number. Values of K c and K g 
are obtained from Figures 5-3, -4, -5 of reference [8]. 



If we consider an order of magnitude approximation, the first term 
in (A-22) is by far the greatest contributor to friction factor for 
small pressure differentials. 

P +P 
12 

The approximation, therefore, that — - — ■ p ra ~ P 1 “ P 2 
reduces (A-22) to the following 

f - [ 2 8 pp m f - < W + ¥ (1+p2) ] ^ <A_23) 

^ m 

London [11] points out that the relationships given here can be 

recombined in such a manner so as to show that f is proportional to the 

3 

porosity cubed and inversely proportional to the compactness (f a p /6). 
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COLBURN 1 FACTOR ; 



Colburn j factor is defined as; 



j = N 0 N„ 

J St Pr 



2/3 



_h_ N 2/3 

Gc Pr 
P 



Substituting (A-17) for G and multiplying by A/A yields 

2/3 



A 

4 _ hA c M 
me A Pr 
P 



but P- = N 

me Tu 
P 



therefore : 



j = N — N 
J Tu A Pr 



Combining equations (A-l) and (A-3) : 



thus : 




BL 



j - N, 



Tu Pr 



2/3 1 * 
L B 



(A-24) 



(A- 25) 



(A-26) 



(A-27) 



It is apparent from equation (A-27) that the Colburn j factor 
is directly proportional to porosity and inversely proportional to 
compactness (jap/B). 



HEAT TRANSFER POWER AND FLOW FRICTION POWER: 



An evaluation of the heat transfer power versus flow friction 
power is of interest in that it is a measure of relative performance. 

The heat transfer power per unit area per degree temperature 
difference is [8] : 



h = 






N 



1 

2/3 D. 



j N 



H 



R 



(A-28) 



pr 

Equation (A-28) with the properties evaluated at standard conditions 
of dry air at 500°F and one atmosphere becomes: 
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h STD = 0.02195 (N R j ) Btu/hr ft 2 °F 

H 

where 

c = 0.2477 Btu/lbm °F 
P 

M = 0.0678 lbra/hr ft 
p = 0.0413 lbm/ft 3 



(A- 29) 



N_ = 0.671 
Pr 

The flow friction power per unit area is [8]: 

3 



E = ^ ( jr} f N 

2g c p ' D h ' 



(A-30) 



When equation (A-30) is evaluated at standard conditions (y and p 
values given above) , the flow friction power per unit area at standard 
conditions is : „ 

E, 



STD 



- 1 ' Ul< 10 ' 7 (t ) 3 '('iSSo') HP/ft: 

n 



For comparison purposes the surface geometries were reduced to a 

-3 

common hydraulic diameter of D = 2 x 10 ft. 

n 
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APPENDIX B 



DESCRIPTION OF EQUIPMENT 

Heat transfer data was obtained utilizing the "single blow 11 
technique which consists briefly of monitoring the fluid temperature 
response at the test matrix exit while subjecting the inlet to the test 
matrix to a step change in fluid temperature* 

Friction factor was determined from pressure drop data obtained 
from static pressure taps located in the test section - one at inlet, 
the other at the exit* 

Necessary equipment to perform such an experiment falls into one 
of the following categories: 



(1) 


Fluid source 




(2) 


Flow metering system 




(3) 


Temperature measuring 


system 


(4) 


Pressure measuring system 


(5) 


Fluid heater system 




(6) 


Matrix holder and test 


section casing 



FLUID SOURCE: 

Air, the working fluid, was provided to the test apparatus by 
connecting the rig to the inlet of a 30HP, multistage, Spencer Turbo- 
Compressor rated for 550cfm operating on a 220 volt a.c. power supply. 

FLOW METERING SYSTEM: (See Figure 1) 

The flow metering system consisted of an ASME standard orifice 
section constructed for d and d/2 pressure taps in a 3.08" inside diameter 
metal tube. [14]. A wide flow range was obtained by utilizing thin 
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plate concentric orifices of throat diameter 2*310, 1*971, 10 232, 

1*540, 0*462, and 0*308 inches respectively* 

Control of the pressure drop across the orifice was maintained by 
a gate valve downstream of the apparatus, ahead of the compressor inlet 
and another gate valve on the compressor inlet * The blast gate on the 
compressor discharge was pre-set so that the unit could not be operated 
beyond the full load rating of the turbo 0 

TEMPERATURE MEASURING SYSTEM: (See Figure 1) 

The majority of the temperatures were measured with iron-cons tantan 
thermocouples. These thermocouples use iron for the positive conductor 
and constantan for the negative conductor., Measuring T^ is a group of 
five thermocouples in series, bound together and insulated from each 
other by teflon tape and inserted in a 1/8 M diameter aluminum tube to 
serve as a radiation shieldo A small aperture was cut in the inserted 
tube facing upstream, T^ and T^ are five wire thermocouple grids, connec- 
ted in series, mounted permanently in the test section casing, T^ is 
a movable five wire series connected thermocouple grid located in the 
matrix holder so that it is adjacent to the downstream side of the 
matrix* The purpose of connecting the thermocouples in series was to 
magnify the emf output so that the instrument sensitivity to small 
temperature changes would be enhanced* With T^ bucked against it 
is possible to determine the temperature uniformity across the matrix* 

Tj versus T^ indicates the difference between upstream and downstream 
temperatures; thus, it is this combination which is used to record the 
rate of change of upstream and downstream temperatures with time, result- 
ing from the step change in fluid temperature upstream. This differential 
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is recorded on a Minneapolis-Honeywell "Brown" Recorder and serves as 
the primary data for heat transfer evaluation of the matrix. The 
recorder has variable chart speeds so that acceptable cooling curves for 
maximum slope determination can be generated. This instrument has a span 
adjustment which permits the span to be varied continuously from the 
narrowest to the widest span desired from 0 to 55 millivolts. For the 
work performed herein, the recorder was precalibrated for a 0 - 3 milli- 
volt scale. This unit also has adjustable supression and damping adjust- 
ment features. The damping adjustment provides a filter network for 
removing a.c. strays, and provides correct damping so that neither under** 
shoot nor overshoot in the response curves is experienced. 

A copper constantan thermocouple was inserted into the duct just 
ahead of the orifice to measure T , orifice temperature. This temperature 
was read in millivolts on a Rubicon Company Portable Precision Potentio- 
meter, converted to degrees Fahrenheit and recorded for each run. 

A copper constantan thermocouple grid was inserted in place of T^ 
to determine the temperature distribution across the cross section during 
initial check out of the apparatus. This thermocouple arrangement was 
part of a special holder lined with balsa wood to provide a smooth, 
continuous flow passage, which also contained an impact tube used for 
checking the velocity profiles (See Figure 33). 

PRESSURE MEASURING SYSTEM; (See Figure 1) 

Pressure taps are located upstream and downstream of the matrix 
and on either side of the orifice. Each tap is connected to an approp- 
riate manometer or draft gage via Imperial Company "poly-flo" tubing. 

Quick closing valves were installed at various positions in the line to 
permit isolation of sections and as a safety feature. 



67 



The following instruments were utilized? 

(1) Meriam Instrument Company, 0 - 50" manometer 

(2) Meriam Instrument Company, -8" to -t-8" manometer 

(3) Ellison Differential Direct Draft Gages, 0 * 6”,0 - 8" 

(4) Ellison Inclined Draft Gages, 0 - 0®5 19 , 0 - 6 ,? ® 

(5) Eo Vernon Hill and Company Type "C" Micromanoraeter , 

0 - 1 o 25 M o 

Any one of these instruments or a combination can be used to measure the 
differential pressures of the orifice or matrix or the required static 
pressures - upstream of matrix and upstream of the orifice® Cross check- 
ing of the various instruments assured reliable operation ® 

FLUID HEATER SYSTEM: (See Figure 1) 

The heater section is comprised of 28 nichrome wire heaters of 
0.0031" diameter 0 The heater system was designed to elevate the air 
temperature 20®F above ambient for a mass rate of flow of 1000 Ib/hour® 

The nichrome material was selected because it has high resistivity, low 
thermal conductivity and specific heat, and has a very small time constant 
so that it permits one to approach the idealized step function® 

The wire heaters are wound two to a bakelite frame, l/32 n between 
each wire, with 50 and 52 wires respectively to each heater on a frame® 
These heaters are connected in parallel, then via a switch to a variable 
voltage bus® The heater frame and its electrical connections are apparent 
in Figures 34 and 35® 

A schematic wiring diagram for one of the two heater complexes is 
shown below® A total of 14 frames are wired in this fashion, thereby 
putting all 28 heaters in parallel, permitting two to be switched on or 
off at a time® 
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The number of heaters in use decreases as flow rate decreases. The 
necessary voltage variations are obtainable by a General Radio Company 
Type W20HM "Variac" Autotrans former, 0 - 280 volt, 8 amp load from a 
240 volt 50 - 60 cycle line. 

MATRIX HOLDER AND TEST SECTION: (See Figures 31, 32, and 35) 

The matrix holder and test section are made of polyethylene plastic. 
The test section casing holds the T^ and T^ thermocouple grids permanently 
mounted, as well as the upstream and downstream static pressure taps 
required for determining friction factor. All parts were machined to 
close tolerance to assure a snug fit for the matrix holder and to guaran- 
tee good alignment of the flow channel through the heater and test sec- 
tions. The flow channel is 3 - 3/16" x 3 - 3/16" and matrices of flow 
lengths up to 3 M may be tested. As mentioned previously under Temperature 
Measuring System, the T^ thermocouple grid is located in the matrix holder. 

Correct positioning of the test matrices in the flow channel is 
ensured by using styrofoam plastic inserts. The inlet cone and stainless 
steel 60 mesh screen employed provided a uniform velocity to the heater 
section. 
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APPENDIX C 



PERFORATED NICKEL GEOMETRIC PROPERTIES 
The material used in the perforated plate matrices was pure nickel, 
electro-deposited sheet of integral structure manufactured by Perforated 
Products, Incorporated** Types 160/40TV, 160/40Q, 50G, 125M and 125P 
were utilized in the matrices - the first two of these having slotted 
openings and the remaining three having round openings 0 A brief des- 
cription of each of the aforementioned types as specified by the manu- 
facturer is set forth in Table C-lo 

Inasmuch as the perforations in all cases were conical, the use of 
the average slot length or average hole diameter was considered a. 
propos , Close examination of each of the different types revealed that 
this average value was somewhat less than the dimensions specified by the 
manufacturer in each caseo In addition to the requirement for determin- 
ing porosity of the plate, a correction for the increased longitudinal 
conduction length due to the perforations as well as a correction to 
the solid conduction cross-sectional area was performedo The following 
idealized patterns for each type used, best illustrate these corrections 0 
160/40TV 
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Cross Hatched Area, = (,0063) (.025) = .0001575 in. 2 
Slotted Area, A gl , in A ch = (.0008) (.0155) = .0000124 in. 2 
Plate Porosity, the fraction of open or slotted area. 



v 



Solidity, s, is 1 - v 



s 




. 0000124 

.0001575 



0.9212 







i 



conduction path length, = 0.0125 + 0.0063 = 0.0188" 




.0063 

.0188 



0.3351 



The solid cross sectional area for heat conduction, A. varies along 

K i 

the path length. The area, Aj^, selected was the cross sectional area 

between the slots normal to the flow direction, i.e., 

A, = (0.0095) (0.0022) = 0.0000209 in. 2 
i 

Just as L = EL^ where L is the total matrix flow length and 
is the flow length between two successive perforations; so does = EL^ 

and A^ = EA^ . In evaluating A^, a conduction area reduction ratio 

j _ cross sectional area of perforations 

is employed and is defined as t — K ; 

cross sectional area 

therefore, A^, is the solid cross sectional area multiplied by 1 - 4 
i.e. , 
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160/400. 



\ = AJ1-C) 



- 1 ' 6058 ( 1 - ( ( P 025 H : o 022 )) - °’ 6102 i ” 2 



= 0.00424 ftc 



r- 



. 02 ? 




A ch = (.0063) (.025) 



,0001575 in. 



A « (.0019) (o0161) * .00003059 in. 
v = .1942 
s = .8058 



,655 



A k - A (1-0 * .002887 ft. 



— =0.3351, same as 160/40TV 
L k 




0.19466 
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. 0079 sin 60 * 
=.0QG84-" 



The equivalent conduction length, , is computed in the case of 

i 2irr k 

round perforations assuming the path to be as shown above, i.e., — — 

where r^ is the radial distance equal to one-half the center to center 



distance. Thus 



_ 2tt (.0079) 
’k 3 2 



.00827" 



and 




.00684 

.00827 



.82709 



is determined as in the case of the slotted perforations. 

_ diam. of hole x plate thickness 

center to center distance x plate thickness 



( . 00366) ( . 0016) 

(.0079) (.0016) * 

Therefore, = A g (1-c) = 1.1679 in. 2 (1-.4633) 

= .6268 in. 2 = .00435 ft. 2 
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125P 






A = (.00684) (.0079) = .00005404 in 2 
ch 

Ahoie = \ ( . 00291) 2 = .00000665 

v = 0.12307 

s = 0.8769 

L, , = .00827" 
ki 



= A g (l-C) = .00640 ft 2 



-- = .82709 
L ki 

C = 0.3684 



50G 




A ch = (.01706) (.0197) - .00033608 in. 2 

Ahoi e - |-(.0138) 2 = .00014957 in. 2 

v = .44504 

s = .55496 

I* = .02063" 

K i 




.82697 



4 - .7005 

A k = A (1-5) s >00243 ft. 2 
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Type Opening Hole Diameter Plate Thickness Center to Center % Open Area 

( inches ) Cinches) Distance, (inches) 
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APPENDIX D 



SAMPLE CALCULATIONS 

The following sample calculations are based on data obtained from 
Run #1 for the modified parallel plate matrix constructed from the 
160/40Q perforated nickel material. The basic geometric parameters are 
illustrated in Appendix A. 

Prior to the reduction of data, several fixed parameters must be 
evaluated. 

DIMENSIONS OF FINNED SHEETS: 3.9375" x 2.00" x .0016" 

(before forming) 

DIMENSIONS OF SPLITTER PLATES: 3.1875" x 2.00" x 0016" 

NUMBER OF FINNED SHEETS: 102 

NUMBER OF SPLITTER PLATES: 103 

WEIGHT OF MATRIX: 311.0979 gms/453.6 gm/lb = 0.68584 lb. 

MATERIAL CONSTANTS: 

k =38.7 Btu/hr - C F - ft 
s 

c = 0.1065 Btu/lb - ®F 
s 

p - 0.321 lb/in 3 - 554.7 lb/Ft 3 

PLANE SURFACE AREA, !*, - ( !0 2) < 2 M 3 . 9 375-P3 .1875? , 

144 

20.1875 ft 2 

HEAT TRANSFER AREA, A, = PLANE SURFACE AREA x SOLIDITY 

A = A*s 

= (20.1875) (.8058) 

A = 16.2667 ft 2 | 

FRONTAL AREA, A fr , = 3.141" x 3.1875"= 10.012 in 2 

A, = 0.06953 ft 2 
fr I 
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It has been shown in Appendix A; however, that this can be alternately 



expressed as 



D = 4 £ 

H 3 

Inasmuch as the perforations have no effect on the hydraulic diameter 
the value of 3 rather than 3 is used. 



_ _ (4) (0.8834) 

H (1742.1) 



0.002028 ft 



Now we are in a position to reduce the data. The following 
presentation outlines the steps required. 

RECORDED DATA: 



d = 2.310" 
o 

AP = 6.69 in H.O 
o 2 

AP = 10.15 in H.O 
m 2 

P = 10.65 in HO 
s 2 

P = 20.80 in H„0 
o l 



T = 69.9 



P = 30.140 in Hg. 
atm 

Chart Speed = 4 in/sec. 
*o= d o /d= -f^r = 0 - 75 



DETERMINATION OF MASS RATE OF FLOW: (See Appendix A) 

P„ 



m = 



where : 



589.81 Fa Y do C 






A P o X x 

o 



C. + AC 
o 



Lvd >] 



from reference [13] 



P = (P - P /13.6) 0.4912; T - t (°F) + 459.7 
3 atm o o o 

Y = 1.0 Fig. 40b of Reference [14] 

F = 1.0 Fig. 30 of Reference [14] 

cl 

For 3„ = 0.75 and d =-3.08", C q = 0.60691 and AC = 0.03839, respectively. 



P 3 " 



30.140 - 



20.80 

13.6 



x 0.4912 = 14.Q537 psi 
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SOLID MATRIX CROSS SECTIONAL AREA, A , 

s * 

A = 102 x 3o9375" x .0016" + 103 x 3, 1875" x .0016" 
s 

= 1.1679 in 2 



A 

s 


= 0.00811 ft 2 


! FLOW AREA, A . - ! 




e 


A 


s A £ - A 


c 


fr s 


A 


= 0.06142 ft 2 


c 





CONDUCTION AREA CORRECTED FOR EFFECTS OF PERFORATIONS, A^ 

A^ = A^ (1-c) where £ is the conduction area reduction 
ratio (See Appendix C) „ 



= 0.00811 (1-0.644) 



\ 



0.00288 ft 



pnnnQTTV MATRIX FLOW VOID VOLUME Vm-Vs _ Vs 

’ P ” MATRIX VOLUME Vtn L ~ Vm 



P = 1 ~ 



0.00135 

0.01159 



0.8834 



rOMPACTNF^ T - TOTAL HEAT TRANSFER AREA 
COMPACTNESS, 6 - matr^ VOLUME 

A 16.2667 



Vm 0.01159 



1403.75 ft 2 /ft 3 



Note that 6 includes the effect of area reduction due to perforations. 
B for an unperforated surface would therefore simply be the plane 

surface area divided by the matrix volume. 

A* 20.1875 ...... . 2 3 

B = — = — - 1742.1 ft /ft 



Vm 0.01159 



THE HYDRAULIC DIAMETER, D„, = 4r u = 

H h 



4 FLOW CROSS SECTIONAL AREA 
WETTED PERIMETER 



T = 69.0 + 459.7 = 528.7 ®R 
o 

• (589.81) (2. 310) 2 

m = y - ^ . . ■■ ■ x 

/l - (0.75) 4 
m = 1605.242 C 






69 x 



14.0537 

528.7 



assume C = 0.6000 
then m = 963.145 Ibm/hr 
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C - C + AC 
o 



10 ' 



N 



R(d o ) 



where 



N 

R(do) 



(4) C12)^ 

(3. 1416)(3. 08)y 

4.961 - 
M 



y is evaluated at temperature t Q by means of the linear approxima- 
tion y = 0.0395 + 0.64167 x 10~ 4 t Q = 0.0395 + 0.64167 x 10~ 4 (69.0) 

= 0.04392 lbm/hr ft 



therefore, 

C = 0.60691 + 0.03839 



(10 4 ) (0.04392) 
(4.961) (963.145) 



C = 0.6104 

This value of C is compared with the assumed value of C and 
successive iterations are performed until the two values are equal* 
The iterative process shows that a value of C = 0*6172 satisfies 
this requirement and thus 

m = 990.756 lbra/hr 



MATRIX REYNOLDS NUMBERS: 



Two values for Reynolds number are required, one for the isothermal 
flow friction factor evaluation and second for the heat transfer 
evaluation* The two values differ only in so far as the difference 
in the absolute viscosities, y^ is evaluated at orifice temperature 
(t ) whereas y„ (heat transfer) is evaluated at average bulk fluid 

O H 

temperature which is assumed to be at t Q + 10°Fo 



N 



Rf 




(4)^(990.756) 
(16.2667) (.04392) 



924.54 



N 



RH 



4Lm 

% 



= N„, — = 924.54 

Rf y H 



(.04392) 

(.04456) 



911.23 
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CONDUCTION PARAMETER: 



A = 



k s\ 



= 0.00282 



Lmc (2/12) (990.756) (.24) 

P 

X, = X7- ; 7— = 0.3351 from Appendix C. 

k L k L k 

therefore 



A k = (.00282) (.3351) = 



0.000944 



MAXIMUM SLOPE: 

The slope of the cooling curve was determined by the method 

prescribed in Appendix A from Figure 36. Slopes were determined 

for each of the curves and the average value was used. 

C ' - ” 



MAXIMUM SLOPE - 



- C 

dll _ _s_ 

d0( " C 

-'max 



x SLOPE x CHART SPEED 



C = W c - (0.68584 lb) (0.1065 Btu/lbm 9 F) = 0.07304 
s s s 

C = me - (990.756 lb /hr) (0.24 Btu/lbm °F)(1 hr/ 3600 sec) = 0.06605. 
P 

SLOPE = 0.14357 (See Figure 36) 



whence , 



MAXIMUM SLOPE = (0.14357) (4) = 



0.63507 



(.06605) 

COLBURN j FACTOR : 

Enter Table 1 or Figure 4 of Howard [5] or Figure 2-A of Howard 

[4] with the values of X, and MAXIMUM SLOPE to determine N_ * 4.05. 

k Tu 

Linear interpolation from this table is sufficient. 

A c 2/3 

j “ n tu r~ n pr from A pp endix A » 

2/3 -4 

N 0.8017 - 0.82353 x 10 t using a linear approximation; 

i K o 



therefore 



j = (A ° 05) (16?2667} ( “ 7960) = 



0.01217 



FANNING FRICTION FACTOR: 



From Appendix A, the friction factor is: 



f ^ 



j~2gep 



6P 
m 



- ( K c -K e) - 
80 



AP 



m 



mean 



(1+P 2 ) 



']? 



14.8048 psi 



P „ = 30.140 in Hg x 0.4912 ft 8 - * = 

atm in Hg 

P = 10.65 in H„0 x 5.204 7P — 5 -~ ■ 55.42 psf = 0.3849 psi 

s 2 in H_0 r 

P = 10.15 in HO x 5.204 . ■ = 52.82 psf = 0.3668 psi 

m 2 in H 2 0 r r 

P = P „ - P = 14.8048 - 0.3849 = 14.4199 psi 

1 atm s 



P„ = P, - AP = 14.4199 - 0.3668 = 14.0531 psi 
A 1 ra 



P +P 

12 



mean 



mean 



= 14.2365 psi 
(144) 



(14. 2365) (144) „ .,.3 

P = p/„ ' ./cn ' T\ - ) ' cV o ' cw ' c ' oo -n = 0.07268 lbm/ft 

mean R(t +459.7) (53. 35) (528. 7) 



G = 



m 



o 

990.756 

0.06142 



= 16130.84 lbm/hr ft 2 = 4.481 lbm/sec ft 2 



fH _ 0.000507 
L 2/12 



0.00304 



K c = 0.48 and K g = -0.33 from Figure 5-3 of reference [8J. 

(1+p 2 ) = 1 + (0.8834) 2 = 1.7804 

- f(2) (32. 2) (0.07268) (52.82) f /R _ .3668 

f L (4.481) 2 (0 * 48 14.2365 



(1.7804) X 



■J 



(0.00304) 



= (12.3126 - 0.150 - 0.04587) (0.00304) 
f = 0.03686 

2/3 

The thermal properties, y and Np^ , were obtained by linear 

interpolation of data from reference [3]. The values of k and c 

s s 

were obtained directly from references [2] and [1] for nickel and 
stainless steel respectively. 
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